Influence of secondary decay on odd-even staggering of fragment cross sections 
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Odd- Even Staggering (OES) appears in many areas of nuclear physics, and is generally associated 
with the pairing term in the nuclear binding energy. To explore this effect, we use the Improved 
Statistical Multifragmentation Model to populate an ensemble of hot primary fragments, which 
are then de-excited using the Weisskopf-Ewing statistical emission formalism. The yields are then 
compared to experimental data. Our results show that, before secondary decay, OES appears only 
in the yields of even mass fragments and not in the yields of odd mass fragments. De-excitation 
of the hot fragments must be taken into account to describe the data, suggesting that the OES in 
fragment yields is a useful criterion for validating or adjusting theoretical de-excitation models. 

PACS numbers: 25.70.Pq,24.60.-k 



I. INTRODUCTION 

Odd-Even Staggering (OES) is a widely observed phe- 
nomenon in nuclear physics. The early observed OES 
in nuclear masses associate the phenomenon with the 
pairing term in the binding energy fl8l . [l9| . Since parti- 
cle production yields, including fission fragments, corre- 
late strongly with the binding energy, measured fragment 
cross sections exhibit OES effects. Sawtooth-shaped 
charge correlations, based on the Z distributions of frag- 
ments produced in different reaction mechanisms includ- 
ing fission and multifragmentation, have been reported 
[T5l - [l7| . Much of the previous work focused on stagger- 
ing as a function of the atomic number, Z, because most 
experimental data was limited to only elemental identifi- 
cation. 

Recently, isotopic identification for heavy reaction 
products can be achieved with state-of-the-art detectors. 
When the charge distributions are subdivided according 
to the neutron excess of the fragments, the staggering 
plots with isotope resolution reveal a more complex struc- 
ture suggesting the de-excitation of the hot fragments 
contribute significantly to the observed OES effect. In 
fragments that are very neutron-rich, the OES effect may 
be reversed, i.e. production of odd mass fragments are 
enhanced, compared to the less n-rich fragments. In this 
work, we investigate the OES effects of the isotopic frag- 
ment distributions obtained in the projectile fragmenta- 
tion of 40 Ca, 48 Ca, 58 Ni, and 64 Ni at 140 MeV/nucleon 
whose experimental analysis has been reported in Rcf. 

Using the grand-canonical version of the Improved Sta- 
tistical Multifragmentation Model (ISMM) [2l| and the 
Weisskopf-Ewing statistical emission to describe the de- 
excitation of the hot primary fragments, we examine the 
effects associated with the latter on the observed stag- 
gering. Our results suggest that the OES effect can be 



very useful in constraining the treatments adopted in the 
description of the de-excitation process. The remainder 
of the manuscript is organized as follows. In Sec. [TT] we 
briefly recall the main features of the ISMM. The results 
are presented and discussed in Sec. IIII1 and the conclu- 
sions are presented in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

In the framework of the grand-canonical approach 
[iH [23j|, the yields Y(A, Z) of a fragment with mass and 
atomic numbers, respectively, A and Z , reads: 



Y(A,Z) 



9A,zV f A^ 



(1) 

where g stands for the spin degeneracy (taken as unit, 
except for the empirical values used for A < 4), At = 

2irh 2 /m n T is the thermal wavelength, m n denotes the 
nucleon mass, fx p (fib) represents the proton (baryon) 
chemical potential, f% m z° th (T) is the Helmholtz free en- 
ergy associated with the fragment plus the pairing term 
of the binding energy B^ l ™ s . The free volume reads 
Vf = x^o-, where Vq is the source's volume at normal 
density and we use x — 2 throughout this work. As in 
Refs. [ill, [23j], the Helmholtz free energy has contribu- 
tions from the fragment's binding energy Ba z, terms 
associated with the Wigner-Seitz correction [24[ to the 

; A1/3 ( 1+ ^i/ 3 , besides those related to 
the internal excitation of the fragment z : 



Coulomb energy C c 
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For clarity, we use the standard SMM parameterization 
of the internal free energy [22j , which is isospin indepen- 
dent, i.e., f\ Z {T) = f\{T). The pairing term is written 
as: 



Y(z)=j2y(Az)/j2 Y (Az) 



(8) 



A,Z 
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(3) 



where N = A — Z denotes the neutron number, and S p 
is the pairing energy. In this work we employ the liquid 
drop mass formula presented in Rcf. 



211, and we refer 



the reader to that work for the numerical values of the 
parameters. Although we have employed a liquid drop 
parametrization of the binding energy in all the calcu- 
lations, we have checked that our conclusions are not 
affected by the use of more precise values of the binding 
energy, such as those given by the procedure described 
in Ref. [HI. 

With these definitions, Y(A, Z) may be rewritten as 
the product of a smooth term multiplied by a rapidly 
oscillating function of Z: 

Y(A, Z) = y smooth (A Z) x e (-D^*[i + (-D N - z ] , (4) 



where 



^smooth (A, Z) = e 



-[fTz(T)-v P Z-v b A}/T 



(5) 



The chemical potentials are determined by simultane- 
ously solving the equations: 



A =^AY(A, Z) 



(0) 



A,Z 



and 



Z Q = J2ZY(A,Z) 



(7) 



A.Z 



where Aq (Zq) denotes the mass (atomic) number of the 
decaying source. 

The de-excitation of the hot primary fragments is 
taken into account through the Weisskopf-Ewing statis- 
tical emission, as described in Refs. [25|,[26|. We consider 
the emission of nuclei up to 20 O. For consistency, the 
same values of the binding energy Ba,z are also used 
in the decay treatment. The same remarks hold for the 
density of excited states, as explained in Ref. [25| . 



III. RESULTS 

The experimental Z distributions obtained in the pro- 
jectile fragmentation of 40 Ca, 48 Ca, 58 NL_and 64 Ni at 140 
MeV/nucleon on 9 Be reported in Ref. [20| are shown in 
Fig. [TJ where the normalized quantity 



is plotted for each case. The data, represented by the full 
circles, increases, on the average, as a function of the frag- 
ment's atomic number Z. A detailed inspection reveals 
that Y(Z) deviates up and down from the local average 
value. As observed before, this staggering is weaker in 
the reactions with more neutron-rich projectiles. 
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FIG. 1. (Color online) Z distribution of fragments observed 
in the projectile fragmentation of several projectiles on a 9 Be 
target. The data are taken from Ref. [2TJ| whereas the calcula- 
tions correspond to the grand-canonical version of the ISMM, 
using the liquid drop formula of Ref. [2l[, at breakup temper- 
ature T = 4.0 MeV. For details, see the text. 

The predictions made by the ISMM, presented in Sec. 
HIl are also displayed in this figure and are represented by 
the full (primary yields) and dashed (final yields) lines. 
There is no attempt to fit the data with the calcula- 
tion. The model is primarily used to give insights into the 
OES. Instead of doing calculations with different values 
of temperature in order to account for the variation of the 
centrality of the collisions, for clarity and simplicity, we 
adopted a single average value of T for all systems. The 
primary distribution was generated using T = 4.0 MeV. 
This choice was based on the criterion that the yields, 
after the de-excitation of the primary fragments, would 
follow the experimental distribution as close as possible. 
Slightly different average temperatures could have been 
used, but this would not affect the conclusions of this 
paper. The model results follow the experimental trend 
of exhibiting weaker staggering effects in the case of the 
neutron rich projectiles. Furthermore, by comparing the 
primary and final distributions, one sees that the devi- 
ations are clearly enhanced by the decay of the excited 
fragments. 

To isolate the local staggering behavior, following Rcf. 
17J, we obtain the average value Y z{Z) by carrying out 



a parabolic fit, considering two points to the left and to 
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the right of Z, besides its own value, except when Z is 
close to or at either ends of the Z range, where the points 
lying beyond the data edge are obviously not used. More 
precisely, although five points are always employed in 
the fit, one shifts the selected region so that it does not 
extend beyond the edges. The ratio: 

Rz{Z) = Y Z {Z)/Y Z {Z) (9) 

is plotted in Fig. [31 for the data from Ref. [2(| shown 
in Fig. [TJ The features of this figure are consistent 
with previous studies. There is a weaker effect for the 
neutron-rich projectile, and the trend generally decreases 
with increasing Z. 




o.o 1 1 

10 20 10 20 30 

atomic number 

FIG. 2. (Color online) Ratio between Y Z (Z) and Y z(Z) for 
the four projectiles in Ref. [2C| . 

In the literature, the staggering effects have been as- 
sociated with the pairing energy, S p . According to Eq. 
[H the staggering effects would be different depending on 
whether the neutron excess, N — Z, is even or odd. In the 
ISMM framework, the only non-smooth contribution to 
the primary yields is the empirical binding energies. This 
is illustrated in Fig. [3] for the 58 Ni system. To explore 
this issue, we subdivide the experimental yields by their 
neutron-excess, and plot even and odd values separately. 

Y N _ Z {Z)=Y{A,Z)/Y,Y{AZ) , (10) 
z 

which gives the yields of an isotope Z with neutron 
excess N — Z. Based on Eq. (j4|), the staggering ef- 
fect should be absent for odd values of N — Z and 
Y N - Z {Z) will be smooth for N - Z odd, since Y(A, Z) 
becomes Smooth (A Z), whereas important staggering ef- 
fects should be expected when N — Z is even due to the 
multiplying factor e^ 1 ^ 3tt. 
The ratio: 

R N -z(Z) = Y N _z{Z)/Y N _z{Z) (11) 



is plotted in the bottom panels of Fig. [3] for different val- 
ues of neutron excess N — Z = 1,3,5 (left panels) and 
N — Z = 0,2,4 (right panels), obtained from the pri- 
mary yields of the ISMM. For succinctness, we focus on 
the 58 Ni projectile, but similar results are also obtained 
with the other projectiles. As expected from Eq. ((3]), the 
primary yields from ISMM clearly show that the OES is 
not obvserved in the case of odd neutron excess values 
whereas the OES is observed for even N — Z values and 
the staggering diminishes as Z increases. Aside from the 
different behavior for odd and even values of N — Z, the 
staggering effects seem to be independent of the N — Z 
values. This behavior is consistent with the empirical 
binding energies, as illustrated in the top panels of Fig. 
[31 where we plot the same type of ratio as above, but 
using the binding energies: 

R(BE) N _z{Z) = BE N - Z {Z)I BE N . Z (Z) (12) 

Fig. E] shows that staggering in the yields of fragments 
with odd neutron excess is not related to the primary 
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FIG. 3. (Color online) Comparison between staggering behav- 
ior in experimental binding energies (top panels) and primary 
yields from the ISMM (bottom panels). Fragments with odd 
neutron excess are shown on the left panels, while fragments 
with even neutron excess are shown on the right. 

Since most of the hot primary fragments have already 
decayed prior to being observed at the detectors, we use 
the implementation of the Weisskopf-Ewing statistical 
emission presented in Refs. [l^, [2(| to estimate the ef- 
fects of the de-excitation process on the fragment yields. 
The corresponding ratios are exhibited in the top panels 
of Fig. [4] One notes that the fragments' de-excitation 
enhances the deviations from the average values so that 
Rn-z now exhibits fluctuations on the same order of 
magnitude as those observed in the experimental data for 
N — Z odd. However, the fluctuations are much larger in 
the case of even values of neutron excess, about twice the 
experimental values. One should note that the scale used 
in the panels corresponding to the model results on the 
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FIG. 4. (Color online) The results obtained using the ISMM 
with Weisskopf-Ewing de-excitation are displayed in the top 
panels, and the experimental results obtained with the pro- 
jectile fragmentation data reported in Ref. [2(| are shown in 
the bottom panels. Odd values of neutron excess are exhib- 
ited in the left panels, whereas even values are displayed in 
the right. For details, see the text. 



right hand side of Fig. @] is different from that adopted 
in the other panels of the figures. Thus, although in the 
framework of the ISMM the de-excitation of the primary 
fragments is absolutely necessary to reproduce the exper- 
imentally observed staggering effects in the case of N — Z 
odd, it leads to fluctuations which are much larger than 
those actually observed in the data for fragments with 
N — Z even. Therefore, it strongly suggests that Rn-z 
is very sensitive to the de-excitation scheme used in the 
model calculations and this observable could be used to 
constrain the treatment for the decay of excited frag- 
ments. 

To examine if the reduced staggering for neutron rich 
nuclei (JV — Z — 5) is also observed in the isotopes with 
odd values of neutron excess, N — Z, we turn to the data 
obtained from the neutron rich projectile, 64 Ni. As shown 
in the top panels of Fig. [SJ R5 (joined by dashed lines) 
exhibits an exact opposite trend of the normal even-odd 
staggering effects; the cross-sections of odd-Z nuclei are 
enhanced relative to those with even-Z. The flip in the 
odd-even staggering for JV — Z = 5 isotopes can also be 
seen in the 58 Ni projectile data plotted in the bottom 
left panel of Fig. SJ In addition, in contrast to the trend 
exhibited in even N — Z isotopes where the staggering ef- 
fects decrease with neutron richness, the magnitude of R§ 
is much larger than R\, for the lighter elements [Z < 14). 
The magnitude of the staggering effects in R 5 decreases 
and becomes similar to R\ at higher Z. 

This behavior cannot be easily explained by the ISMM 
since the primary yields lead to Rn~z = 1 for N — Z 
odd, as shown in the bottom left panel of Fig. [3J In 
the framework of this model, it can only be explained by 
intricate correlations associated with the dc-excitation of 



the primary fragments, as one sees in the top right panel 
of Fig. [5j which shows Rn-z computed with the final 
yields. However, the reversed OES takes place only at 
small Z values, while it is observed over the whole Z 
region. 

One striking feature observed in the experimental data 
is the amplification of the staggering in the distribution 
associated with neutron deficient fragments (JV — Z = 
— 1) compared with the ratios obtained with neutron rich 
ones (JV — Z = 1). The corresponding experimental ob- 
servations are shown in the bottom left panel of Fig. [5J 
for the 58 Ni projectile. This amplification has no expla- 
nation in the framework of the ISMM since Eq. (|4| pre- 
dicts that the primary yields would be strictly smooth 
for N — Z = ±1. To check whether the de-excitation 
treatment employed in this work could account for the 
observed enhancement, the ratios obtained with the final 
ISMM yields are displayed in the bottom right panel of 
Fig. [5j This enhancement is not seen in the model calcu- 
lations. Clearly, more detailed investigations are needed 
to understand this feature. 

The theoretical model employed in this work could 
not quantitatively explain many features observed exper- 
imentally. Further theoretical investigation might lead to 
an improvement in the understanding of the de-excitation 
process. For example, this model currently assumes that 
the pairing term in the binding energy has no tempera- 
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FIG. 5. (Color online) Comparison of the ratios obtained with 
N — Z = —1 and JV - Z = 1 for the 58 Ni projectile (lower 
panels) and with JV — Z — 1 and JV — Z — 5 in the case of the 
64 Ni projectile (upper panels). The experimental ratios are 
calculated using the yields of Ref. [23|. For details, see the 
text. 



IV. CONCLUDING REMARKS 

We have examined the effects of the de-excitation on 
the staggering observed in the yields of the fragments 
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produced in projectile fragmentation. The model qual- 
itatively explains this feature through odd-even effects 
associated with the pairing term of the nuclear bind- 
ing energy. However, we have found that the ISMM is 
not able to reproduce the experimentally observed odd- 
even staggering effects in the case of fragments with odd 
N — Z, if the effects associated with the de-excitation of 
the hot primary fragments arc not taken into account. 
Even though the sequential decay fixes this shortcom- 
ing of the model, it leads to fluctuations which are much 
larger than observations in the case of fragments with 
N — Z even. Further examination of the staggering effect 
exposes additional deficiencies of the models employed 
in this work, as the enhancement of the ratio obtained 
with neutron deficient (N — Z = —1) fragments com- 
pared to that calculated with neutron rich (N — Z = 1 ) 
nuclei cannot be reproduced theoretically. Moreover, the 



phase shift in the staggering behavior for TV — Z = 1 and 
N — Z = 5 leads to a flip at large Z. This is only partially 
reproduced by the model, which predicts a flip at small 
Z. These findings reveal that the observed staggering 
is quite sensitive to the treatment employed in describ- 
ing the decay of the primordial hot fragments and can 
therefore be used to constrain these treatments. 
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